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Abstract. We study the dependence structure of market states by estimating 
empirical pairwise copulas of daily stock returns. We consider both original 
returns, which exhibit time-varying trends and volatilities, as well as locally 
normalized ones, where the non-stationarity has been removed. The empirical 
pairwise copula for each state is compared with a bivariate K-copula. This copula 
arises from a recently introduced random matrix model, in which non-stationary 
correlations between returns are modeled by an ensemble of random matrices. 
The comparison reveals overall good agreement between empirical and analytical 
copulas, especially for locally normalized returns. Still, there are some deviations 
in the tails. Furthermore, we find an asymmetry in the dependence structure 
of market states. The empirical pairwise copulas exhibit a stronger lower tail 
dependence, particularly in times of crisis. 
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1. Introduction 

The concept of copulas was introduced by Sklar in 1959 mi to study the linkage 
between multivariate distribution functions and their univariate marginals. Since then, 
copulas have gained growing importance as a tool for modeling statistical dependence 
of random variables in many fields. In finance, the usage of copulas is relatively new 
but it has already found application in risk management, see e.g., 012 010EH3, 
derivative pricing, see e.g., mu Hina Banana, and portfolio optimization, see e.g., 
Gann nans- For an overview of the literature on applications of copulas in finance 
the reader is referred to mm- Copulas allow to separate the dependence structure 
of random variables from their marginal distributions. This is sometimes useful in 
statistical applications as the dependence structure and the marginal distributions can 
be modeled separately and joined together resulting in new multivariate distributions 
with different behavior. For a discussion on difficulties in the application of copulas 
the reader is referred to [2211231 122 • Here, we simply view the copulas as providing a 
standardized way for the study of statistical dependences. The marginal distributions 
are mapped to the uniform distribution; the statistical dependence is considered in 
terms of the marginal cumulative distribution functions. 

Recently, we identified market states as clusters of similar correlation matrices 
and studied their corresponding correlation structures [25] . The correlation structure, 
however, does not capture the full statistical dependence between financial return time 
series. Here, we choose a copula approach to study the dependence structure of market 
states |2f'i| . To this end, we estimate empirical copulas for many stock return pairs and 
average over all of them to obtain an empirical pairwise copula for each market state. 
We stress that the identification of market states relies on the correlation matrices, the 
copulas are used only for analyzing the states and not defining them. To estimate the 
empirical copulas we use both original and locally normalized returns. The original 
return time series exhibit time-varying trends and volatilities E2EH1I22]- These have 
to be taken into account when transforming the marginal distributions to uniform ones. 
To this end, we apply the method of local normalization [30], which leads to stationary 
time series while preserving the correlations between them. The resulting empirical 
copulas provide different information. The copulas for the original returns describe 
the dependence structure on a global scale, i.e., for the full time horizon, whereas the 
copulas for the locally normalized returns describe the dependence structure on a local 
scale. 

The empirical pairwise copulas for each market state are compared with a 
bivariate K-copula, which arises from a random matrix approach introduced in EH [sa¬ 
lt models the non-stationarity of true correlations by an ensemble of random matrices. 
The model yields a multivariate return distribution in terms of a modified Bessel 
function of the second kind, a so-called K-distribution. In [25] the K-distribution was 
found to provide a good description of the heavy-tailed empirical return distributions 
for each market state. Here, we aim to arrive at a consistent description within 
the random matrix model studying the agreement between K-copula and empirical 
dependence structure for each market state. In addition, our study provides further 
evidence for asymmetric dependencies between financial returns [33] 1341 135| . We find 
an asymmetry in the tail dependence of empirical pairwise copulas which we study in 
more detail. 

The paper is organized as follows. In section [2] we review the basic concepts 
of copulas, stating the main result in the copula theory, the Sklar’s theorem, which 
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we use to derive a K-copula. In section [3] we present the data set and recapitulate 
the identification of market states for the Nasdaq Composite market in the period 
1992 — 2013 as done in [25] . In section [4] we study the empirical copula densities for 
each market state and compare them with the K-copula. We conclude our findings in 
section [5] 


2. Copula 


We begin with a short introduction to the concept of copulas in |2.1| For more details 
with an emphasis on the statistical and mathematical foundations of copulas the reader 
is referred to the textbooks of Joe [36] and Nelsen [37]. In 2.2 we present the K-copula 
which plays a central role in this study. We restrict ourselves to the bivariate case, 
since later on we study empirical pairwise copulas. 


2.1. Basic Concepts 

Consider two random variables X and Y. The joint distribution of X and Y contains 
all the statistical information about them. It can be expressed either in terms of the 
joint probability density function (pdf) fx,y{x,y) or in terms of the joint cumulative 
distribution function (cdf) Fx t y{x,y), where 

x y 

Fx,y{.x,y) = J dx' J d y' f x ,y(.x',y') . (1) 


From the joint pdf fx.y( x - y) one can extract the individual distributions of X and Y 
as follows 

OO 

fx(x) = J dy fx,y{x,y) , (2) 

— OO 

and analogously for Y. The densities fx{x) and fy{y), called marginal probability 
density functions, and the corresponding marginal cumulative distribution functions 
Fx( x) and Fy(y) describe the individual statistical behavior of the random variables. 

When dealing with correlated random variables, one is interested in their 
statistical dependence. The Pearson correlation coefficient is commonly used as a 
measure of dependence. It is defined as 

r Cov(X, Y) 

oj,y = - , foj 

(Txcry 

where Cov(X, Y ) is the covariance of both random variables and <Jx and ay are the 
respective standard deviations. However, the correlation coefficient only measures the 
linear dependence between the random variables. 

Copulas provide a natural way to study the statistical dependence of random 
variables. The transformation 


Ui=Fi(i) i = X,Y (4) 

leads to new random variables with uniform distributions on the unit interval, called 
the rank of X and Y, respectively. Their joint distribution is called a copula. It 
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describes the dependence structure of the random variables A' and Y separated from 
their marginal distributions. 

A central result in the copula theory is Sklar’s theorem which enables us to 
separate any multivariate distribution function into two components: the marginal 
distributions of each random variable and their copula 


Fx,y{x,v) = Cop XY (Fx(x),F Y (y)) . 


(5) 


If the marginal distribution functions are continuous, the copula that satisfies 
equation ([5| is given by 

Cop X Y {u,v) = F X:Y (Fx 1 (u), Fy 1 ^)) , (6) 


where F y 1 and F Y 1 represent the inverse cumulative distribution functions, the so 
called quantile functions. This equation allows to extract the dependence structure 
directly from the joint distribution function. From the copula © one can compute 
the copula density as follows 


COP X ,y( u , v ) = °Px,y( u i v ) ■ 


(7) 


2.2. K- copula 


The K-copula arises from a random matrix model introduced in m to model time- 
varying correlations between financial time series [35} |35j. It was first used in [26], 
where it was found to describe the empirical dependencies in financial data much 
better than a Gaussian copula. 

Consider a market consisting of K stocks. At each time t(t = 1,..., T) we assume 
that the return vector r(t ) = (r± (t ),..., r^(t)) is drawn from a multivariate normal 
distribution with a covariance matrix £ t 


g(r\z t ) 


1 

, — exp 

Vdet 2ir £ t 



( 8 ) 


where we suppress the argument t of r to simplify our notation. We now model the 
time-dependent covariance matrix £ t by a Wisliart random matrix AA\ where the 
K x N model matrix A is drawn from a Gaussian distribution with the pdf 

w(A\E,N) = exp trA^-Mj . (9) 

Here, £ represents the average covariance matrix estimated over the sample of all 
r(t),t = 1 ,...,T. Hence, the time-dependent covariance matrices are modeled by 
an ensemble of Wishart matrices AAA which fluctuate around the sample average £. 
Averaging the multivariate normal distribution (8) with the random covariance matrix 
AA* over the Wishart ensemble leads to a K-distribution for the multivariate returns 


($)(r|£, AT) = J d[A] w (A|£,iV) g{r\AA*) 

OO 

1 




(27r)*'T(Af/2)v / det£ 


dz z 2 e 


7rN 


-K 


exp 


r t £' 1 r 

4z J (10) 


V?- N VN K (ViVrt£- 1 r) 

T(IV/2)v'det(27r£) J Nr ^- h 


K-N > 
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where AT„ is the modified Bessel function of the second kind of order v = (K — N)/2. 
It depends only on the average covariance matrix E = aCcr estimated over the whole 
sample, where C is the average correlation matrix and a = diag (ay, ..., ax) the 
diagonal matrix of the standard deviations, and a free parameter TV, which governs 
the variance of the Wishart ensemble 


a a i" i \ E A' / I E k k E 11 

var([AA T ] fci ) = - u 


TV 


( 11 ) 


where E^ is the kl- th element of the average covariance matrix E. Thus, it 
characterizes the strength of fluctuations around the average covariance matrix in the 
considered sample. The larger TV, the smaller the fluctuations around E, eventually 
vanishing in the limit TV —> oo. 

The K-copula is the dependence structure which arises for the K-distribution (10). 
For the bivariate case K = 2, the pdf of the vector r = (n, r 2 ) reads 


fc,N(ri,r 2 ) = (ff)(r|E,TV) = 


1 


T(TV/2) 


■ z N / 2 ~ 1 e~ z TV 

d 2 — , -— exp 


Vl — ' 


Here, we used the covariance matrix 

/ ^2 


E = 


erf erier 2 c^ 
Vcri cr 2 C 






47T2; 


1 c 

c 1 


TV r\ — 2cr\T2 + r\ 


4 z 


1 — c 2 


( 12 ) 


(13) 


where c denotes the average correlation coefficient, estimated over the whole sample. 
We chose the standard deviations one, o\ = <r 2 = 1, since the copula is independent 
of the marginal distributions. Then, the marginal distribution densities are identical, 
fi{n) = / 2 (r 2 ), where 

OO OO . - 

/,(n) = J dr 2 fMn.n) = J iz zf^z-yj A exp (-Er?) . (14) 


According to equation |6]) the bivariate K-copula is given by 
Cop CiN (u,v) = K c , iv (T ? r 1 (M),K 2 _1 ('(;)) , 


(15) 


where c and TV are the parameters of the copula, F c ?? denotes the cumulative 
distribution function of the bivariate distribution (12), and F^ 1 the inverse 
distribution function of the marginal cdf given by 


' l r 2 

F c , N (r i,r 2 )= J d£ j dC f c , N (tC) 


— OO 

ri 


— oo 
ri 


= dU dC 


dz 


z N/2-l e -z N 


T(TV/2) 4irz 


exp 


TV £ 2 - 2c£C + C 5 


4z 


1 — c 2 


= dU dC 


TV 


N(e-2c^+C 2 

l-c 2 


7T r(TV/2)v / 2 JV vTiCc2 


AT 2-i> 


/TV(e 2 -2c^C + C 2 ) 

1 — c 2 


(16) 
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and 


' 1 ' 1 I - 

«(n) = / d f MO = / df^ /d* z N/2 ~ l e~ z \j A exp (-^) 


—oo 
7*1 


= d€ 




Fk r(7V/2) v / 2’ 


JV-1 


/Cl-, 


(v^?) ■ 


(17) 


The K-copula density can be obtained from the K-copula (15) by differentiation 

d 2 Cop c N (u, v) _ f c N [F^ 1 {u)^ F^ 1 (v)) 


cop C , N { U , V ) = 


dudv 


fi(F^(u))f 2 (F^(v)) 


(18) 


It depends only on the average correlation c and the free parameter N, which 
characterizes the fluctuations around E. Figure [l] shows the K-copula density for 
different parameter values. The stronger the average correlation c and the lower the 
parameter N, the higher is the probability for extreme co-movements. 

Furthermore, we note that the K-copula is a symmetric copula. It is based on 
the elliptical distribution (101 and thus it belongs to the class of elliptical copulas. 


3. Identification of market states 


We now present the data set and identify market states as clusters of correlation 
matrices with similar correlation structures, as done in [ 251 . These market states will 
be the object of our empirical study in the next section. 

We consider K = 258 stocks of the Nasdaq Composite index traded in the 22-year 
period from January 1992 to December 2013 [40] which corresponds to 5542 trading 
days. For each stock k we calculate the return time series 


r k (t) = - 


S k (t + At)-S k (t ) 


S k (t) 


k = 1,..., K 


(19) 


where S k (t) is the price of the fc-th stock at time t and At is the return interval which 
we chose to be one trading day. 

Empirical return time series exhibit time-varying drift and volatilities. The 
correlation coefficient © averages over these time-dependent parameters which results 
in an estimation error of the correlations. In order to eliminate this kind of error we 
employ the method of local normalization i50j . For each return time series k we 
subtract the local mean and divide by the local standard deviation 


r{t) - (r{t)) n 
~ y/(r 2 (t)) n -(r(t))l 


( 20 ) 


where (... ) n denotes the local average which runs over the n most recent sampling 
points. For daily data we use n = 13 as discussed in [30] . The local normalization 
removes the local trends and variable volatilities while preserving the correlations 
between the time series. 

Using the locally normalized daily returns we now obtain a set of 131 correlation 
matrices estimated on disjoint two-month intervals of the 22-year observation period. 
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Figure 1: K-copula densities cop cN (u,v) for different parameter values. Top: N = 5, 
(left) c = 0 and (right) c = 0.5 , bottom: c = 0.2, (left) N = 3 and (right) N = 30. 


To identify market states we perform a clustering analysis based on the partitioning 
around medoids algorithm eh, where the number of clusters is estimated via the gap 
statistic [123 . The clustering analysis separates the set of 131 correlation matrices 
into six groups based on the similarity of their correlation structures. Each group 
is associated with a market state. We point out that the identification of market 
states is performed ex-post, the clustering algorithm has the correlation matrices of 
all times. Figure [2] shows the time evolution of the six market states. We observe 
that the market switches back and forth between states. Sometimes it remains in a 
state for a long time, sometimes it jumps briefly to another state and returns back or 
evolves further. On longer time scales the market evolves towards new states, whereas 
previous states die out. 

In the next section we will study the statistical dependence for each market state. 
To obtain the return time series for each state, we proceed as follows: We take the 
complete return time series, r(t) or r(t), and divide it into a sequence of disjoint two 
month intervals. We merge all intervals belonging to a given state according to the 
cluster analysis. We note that the return time series for the six market states differ 
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Figure 2: Time evolution of the market in the observation period 1992 — 2013. Each 
point represents a correlation matrix measured over a two-month time window. 


in length. 

4. Empirical results 

We present the empirical pairwise copula densities for each market states in |4.1| and 
compare them with the K-copula densities in |4.2| In |4.3| we study the asymmetry of 
the tail dependence of the empirical copula densities in more detail. 

4-1- Empirical pairwise copulas for each market state 

To estimate the empirical pairwise copula of two return time series r k {t ) and n(t), we 
first have to transform them into uniformly distributed time series. To achieve this, 
we employ the empirical distribution function 

T 

u k (t) = F k (r k (t )) = 1 { ? ’fc( r ) ^ r fc0)} - > ( 21 ) 

T—l 

where 1 is the indicator function, T denotes the length of the time series and the 
factor 1/2 ensures that the values of the transformed time series u k (t) lie in the 
interval (0,1). The empirical copula density of the time series r k (t) and r;(t) is then 
the two-dimensional histogram of the transformed time series u = u k (t) and v = ui (t). 

An accurate estimation of the copula density requires a large amount of data. 
Thus, for each state we compute the copula densities of all K{K — 1) /2 stock pairs as 
two-dimensional histograms of the transformed time series and then average over all 
pairs 

2 K ~ 1 K 

cop {l \u,v) = K(R _ 1) H cop k}( u > v ) > * = !,••■>6, (22) 

' ' fc=l 2=fc+1 

where the upper index i denotes the state number. For the bin size of the histograms 
we choose A u = Av = 0.05. We note that as the length of the time series for each 
market state is rather short, we are not able to study the full copulas for each stock 
pair k, l separately. Hence, we cannot make any direct statements about similarity and 
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dispersion regarding the full dependence structure. Our results only yield statements 
about the empirical dependence structure on average. 

Figure [3] shows the empirical pairwise copula densities for the original returns 
(191 for each of the six market states. We observe a variation of the dependence 
structure from state to state, particularly visible in the tails. In state 1, which covers 
the period from 1992 to roughly 2000, we find a rather flat copula density, indicating 
low dependence between return pairs. In state 2 we observe deviations from the flat 
copula density particularly in the tails, which become more and more pronounced in 
state 3 and 4. State 5, first appearing during the financial crisis in 2008, exhibits the 
strongest dependence. The dependence decreases again in state 6. 

Fig ure [4| shows the empirical pairwise copula densities for the locally normalized 
returns (20). The dependence structures of the six states are mostly preserved after 
performing local normalization. Deviations are observed in the lower-left and the 
upper-right corners where the copula densities for the original returns exhibit higher 
peaks. The reason for that is the time-varying volatility. During periods of high 
volatility stocks tend to have large returns which contribute to the corners of the 
copula density. On the other hand, in periods of high volatility the correlations 
between stocks become stronger. This leads to higher peaks in the corners of the 
dependence structure for the original returns. 

It is important to note that the copulas for the original and the locally normalized 
returns contain different statements. The copulas for the original returns describe the 
dependence structure for the full time horizon. On the other hand, the copulas for the 
locally normalized returns provide information about the statistical dependence on a 
local scale. 

Furthermore, we observe that the empirical copula densities are asymmetric with 
respect to opposite corners. We find stronger dependence in the lower tail than in the 
upper one, that is, the dependence between large negative returns is stronger than the 
dependence between large positive ones. This asymmetry is an important feature of 
empirical copula densities and thus we will discuss it in more details in section |4.3[ 


4-2. Comparison with the K-copula 

We now compare the empirical pairwise copula densities for each market state with the 
K-copula densities. Again, we consider both the original and the locally normalized 
returns. The K-copula density is obtained in the following way: We calculate the 
K-copula according to equation (15) where the integrals are performed numerically. 
The K-copula density for each bin of size Au = Av = 0.05 is then estimated by 


C0 Pc,n( u i v ) 


Cop e , N (u, v) - Cop gJV (u, v - Av) - Cop s N (u - Au, v ) 

+ Cop gAr (it — Au, v - Av) . 


(23) 


For comparison we compute the difference between empirical and analytical copula 
density for each state 

cop w (u,u) - cop^Oiqt;) , * = 1,..., 6 , (24) 

where c is the average correlation coefficient of all K{I\ — l)/2 stock pairs for the 
considered state. The free parameter N is estimated by a fit which minimizes the mean 
squared difference between empirical and analytical copula density. The parameter 
values for each market state are summarized in table [T] The differences between the 
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Figure 3: Empirical pairwise copula density cop^(u,i;) for the original returns. 
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state 1 state 2 



0 1 2 3 4 5 6 

Figure 4: Empirical pairwise copula density cop*'*) (u, v) for the locally normalized 
returns. 
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returns 


state 1 

state 2 

state 3 

state 4 

state 5 

state 6 

original 

c 

0.046 

0.13 

0.17 

0.25 

0.42 

0.22 


N 

41.7 

11.7 

8.4 

5.6 

2.8 

10.0 


loc. normalized 

c 

0.048 

0.13 

0.19 

0.28 

0.43 

0.25 


N 

70.7 

28.6 

29.8 

15.8 

7.4 

20.4 


Table 1: Parameters of the K-copula density for the original and the locally normalized 
returns. 


returns 

state 1 

state 2 

state 3 

state 4 

state 5 

state 6 

original 

0.11 

0.41 

1.22 

2.47 

5.53 

2.84 

loc. normalized 

0.047 

0.38 

0.70 

1.43 

2.94 

1.50 


Table 2: Mean squared differences between empirical and K-copula densities. 


empirical copula density and the K-copula density for each state are presented in 
figure [5] for the original and in figure [6] for the locally normalized returns. Overall, 
we find a good agreement between empirical and analytical copula densities. The 
K-copula seems to capture the dependence structure of the first three states very well. 
Small deviations from the K-copula density are observed for state 4 and 6. Only the 
dependence structure of state 5 cannot be captured by the K-copula. For the locally 
normalized returns we find a better agreement, which is reflected in the smaller mean 
squared differences, see table [2] Deviations are observed mainly in the corners of 
the copula densities. The empirical copula densities exhibit stronger dependence in 
the tails, i.e., the probability for extreme co-movements is underestimated by the 
K-copula. 

It is important to note that the K-copula captures the empirical dependence 
structure rather well with only one free parameter. However, due to its symmetric 
nature it cannot account for the asymmetry observed in the data. The skewed 
Student’s t-copula is an alternative proposed by [43] which was found to account 
for asymmetric dependencies in financial data |441 I45j . It captures the empirical 
dependence structure of the original returns better than the K-copula due to 
the presence of an additional parameter which accounts for the asymmetry [261 . 
Nevertheless, here we confine ourselves to the comparison of the empirical copulas 
with the K-copula as we aim to arrive at a consistent picture within the random 
matrix model. We discuss the asymmetry in more detail in the following. 

4-3. Asymmetry of the tail dependence 

Asymmetric dependence between returns has been reported by several authors, see 
e.g., ^2] E31 S3- Our study provides further evidence revealing a stronger lower tail 
dependence in the empirical copula densities of market states. We now take a closer 
look at this asymmetry. To this end, we estimate the tail dependence in the four 
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state 1 state2 



-10 12 

Figure 5: Difference between the empirical copula density and the K-copula density 
cop ( l \u,v) — cop^ N (u,v) for the original returns. 
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state 1 state2 



-10 12 

Figure 6: Difference between the empirical copula density and the K-copula density 
cop ( l \u,v) — cop^ N (u,v) for the locally normalized returns. 
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corners for all K(K — l)/2 empirical pairwise copulas according to 


0.2 0.2 



o o 
1 0.2 


ULj^j = J du J dv cop kl(u,v) , 

0.8 o 

l l 

UUg = J du J dv cop kj(u,v) , 

0.8 0.8 
0.2 1 

LU M = J du J dv copjj$(u,v) , 
0 0.8 


(25) 


where the upper index i = 1,..., 6 denotes the state number and the lower indices 
represent a stock pair k,l. Here, LL and UU refer to the lower-lower and upper- 
upper corners, respectively, and represent the positive tail dependence, whereas UL 
and LU refer to the upper-lower and lower-upper corners, respectively, and represent 
the negative tail dependence. The asymmetry in the tail dependence can now be 
quantified by the differences 


d3= UU S-LLg 

0£! = lu 8-ul S 


* = 1,... ,6 


(26) 


where a k \ captures the asymmetry of the positive and /3 k ] the asymmetry of the 
negative tail dependence for each stock pair k,l. Figure [7] shows the histograms of the 
asymmetry values exemplarily for state 5. For the returns we find a negative offset for 
the values of a k { , whereas the values of p k j are centered around zero. This indicates, 
on average, an asymmetry in the positive tail dependence, i.e., simultaneous large 
negative returns are more likely to occur than simultaneous large positive returns. On 
the other hand, we do not find such an asymmetry in the negative tail dependence. 
For the locally normalized returns, we find much weaker asymmetry in the positive tail 
dependence and once again no asymmetry for the negative tail dependence. We note 
that the means of the asymmetry values are the relevant quantities, the asymmetry 

values for each pair are distributed around the mean due to statistical fluctuations. 

(5) (5) 

The standard deviation for a k [ is 0.01 and for j3 k / 0.007. Indeed, the asymmetry 
effect is very small. Still, it is clearly visible, see figure [3| 

In the following, we study the asymmetry values for each market state. Figure [8] 

shows the mean asymmetry values a W for each market state, obtained by averaging 
(i) 

over all a k \ for a given state. On a local scale the asymmetry in the positive tail 
dependence is much weaker. Only state 5 still exhibits a certain amount of asymmetry. 
On the other hand, the asymmetry in the negative tail dependence is negligibly small 
for both original and locally normalized returns. 
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Figure 7: Histograms of the asymmetry values for all stock pairs k,l, exemplarily for 
state 5. Left: asymmetry for the positive tail dependence a k ,, right: asymmetry for 

the negative tail dependence: p k I . Top: for returns, bottom: for locally normalized 
returns. The dashed red lines represent the corresponding mean values. 


5. Conclusion 

We studied the dependence structure of market states by means of a copula approach. 
To this end, we estimated the empirical pairwise copulas for each state and compared 
them with the bivariate K-copula. The bivariate K-copula arises from a random 
matrix model where the non-stationarity of correlations is taken into account by an 
ensemble of random matrices. It is a symmetric, elliptical copula, which depends 
on two parameters: the average correlation coefficient, estimated over the considered 
sample of returns, and a free parameter which characterizes the fluctuations around 
the average correlation in this sample. We estimated the empirical pairwise copulas 
for both original and locally normalized returns. The local normalization removes the 
time-varying trends and volatilities while preserving the dependence structure. The 
corresponding copula describes the dependence structure on a local scale, whereas the 
copula of the original returns provides information about the dependence structure 
on a global scale, i.e., for the full time horizon. Overall, the K-copula captures the 
empirical dependence structure of market states. We found a good agreement, in 
particular for the copulas estimated on a local scale. Thus, we obtain a consistent 
description within the random matrix model: The K-distribution describes the heavy¬ 
tailed return distributions while the K-copula captures the corresponding dependence 
structure. However, we found an asymmetry in the positive tail dependence, i.e., 
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Figure 8: The mean asymmetry values for each state. Left: asymmetry for the positive 
tail dependence aW, right: asymmetry for the negative tail dependence: Top: 

for returns, bottom: for locally normalized returns. 


a stronger lower tail dependence, indicating a larger probability for simultaneous 
extreme negative returns. This asymmetry cannot be captured by our model. It is 
more pronounced on a global scale. On a local scale we find a much weaker asymmetry. 
However, in times of crisis the asymmetry is still clearly present. 
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